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ABSTRACT 

The Canada-France High-z Quasar Survey (CFHQS) is an optical survey designed to locate quasars during 
the epoch of reionization. In this paper we present the discovery of the first four CFHQS quasars at redshift 
greater than 6, including the most distant known quasar, CFHQS J2329-0301 at z = 6.43. We describe the 
observational method used to identify the quasars and present optical, infrared, and millimeter photometry and 
optical and near-infrared spectroscopy. We investigate the dust properties of these quasars finding an unusual 
dust extinction curve for one quasar and a high far-infrared luminosity due to dust emission for another. The 
mean millimeter continuum flux for CFHQS quasars is substantially lower than that for SDSS quasars at the 
same redshift, likely due to a correlation with quasar UV luminosity. For two quasars with sufficiently high 
signal-to-noise optical spectra, we use the spectra to investigate the ionization state of hydrogen at z > 5. For 
CFHQS J 1 509- 1 749 at z = 6 . 1 2, we find significant evolution (beyond a simple extrapolation of lower redshift 
data) in the Gunn-Peterson optical depth at z > 5.4. The line-of-sight to this quasar has one of the highest 
known optical depths at z ~ 5.8. An analysis of the sizes of the highly-ionized near-zones in the spectra of two 
quasars at z = 6.12 and z = 6 A3 suggest the IGM surrounding these quasars was substantially ionized before 
these quasars turned on. Together, these observations point towards an extended reionization process, but we 
caution that cosmic variance is still a major limitation in z > 6 quasar observations. 

Subject headings: cosmology: observations — quasars: general — quasars: emission lines — quasars: absorption 
lines — intergalactic medium 



1. INTRODUCTION 

The earliest observed epoch of the universe is 380 000 years 
after the Big Bang when protons and electrons combined into 
neutral hydrogen. The radiation from this epoch, observed as 
the Cosmic Microwave Background (CMB), shows that the 
universe was homogeneous with fractional density fluctua- 
tions of only ~ 10 -4 . For the following few hundred million 
years there were no stars or galaxies. This period is often re- 
ferred to as the Dark Ages because of the absence of sources 
of light (Rees 1999). 

Eventually these small density fluctuations in the dark mat- 
ter distribution grew into dark matter halos where gas could 
cool and form galaxies. As stars and black holes formed 
within galaxies they emitted copious amounts of ultraviolet 
(UV) radiation. Many of these photons escaped their host 
galaxies and ionized hydrogen in the intergalactic medium 
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(IGM). When sufficient sources of ionizing photons had 
formed, they provided an ionizing background radiation field 
high enough to ionize all the diffuse neutral hydrogen in the 
IGM. The details of reionization are sensitive to the luminos- 
ity function of ionizing sources, the fraction of ionizing pho- 
tons that escape from them, the ionizing source spectra and 
the dumpiness of the IGM (Miralda-Escude et al. 2000). 

Quasars are useful probes of reionization because they can 
be detected at very high redshift and have strong intrinsic UV 
radiation that is absorbed at the resonant Lyman lines by neu- 
tral hydrogen. The presence of significant flux in the Lya 
forest of quasars at z < 6 indicates that the IGM is almost 
completely reionized by this epoch (Fan et al. 2000; Fan et 
al. 2001; Songaila 2004). The discovery of Gunn-Petersen 
troughs in the spectra of quasars at 6 < z < 6.42 have been 
interpreted as evidence that we are approaching the epoch of 
reionization (Becker et al. 2001; White et al. 2003; Fan et 
al. 2006a). Fan et al. (2006b) analyzed a sample of 19 SDSS 
quasars at z > 5.7 and found that the mean Gunn-Peterson 
optical depth and its variance increase sharply away from the 
low redshift extrapolation at z > 5.7, consistent with the over- 
lap phase of reionization occurring at z ~ 6.1 (Gnedin & Fan 
2006). However, it has also been argued that these data are 
compatible with much earlier reionization (Songaila 2004; 
Becker et al. 2007). 

The sizes of highly-ionized near-zones around quasars are 
sensitive to the neutral fraction of hydrogen before the quasar 
turned on (Cen & Haiman 2000). However, the interpretation 
of these results is complex and various authors have come to 
a wide range of conclusions on the basis of the same SDSS 
quasar spectra (Wyithe & Loeb 2004; Wyithe, et al. 2005; Yu 
& Lu 2005; Fan et al. 2006b; Bolton & Haehnelt 2007a,b; 
Mesinger & Haiman 2007; Alvarez & Abel 2007). Despite 
these uncertainties, this could be a powerful probe at high neu- 
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tral fractions where traditional Gunn-Peterson measurements 
are insensitive. Becker et al. (2006) studied O I absorption up 
to z = 6.42 and did not find an O I forest suggesting that these 
lines of sight are either largely unenriched or significantly ion- 
ized. 

Simulations predict that reionization is an inhomogeneous 
process due to the existence of discrete, clustered ionizing 
sources and an inhomogeneous IGM (e.g. Barkana & Loeb 
2001). Cosmic variance in transmission due to large-scale 
density inhomogeneities is extremely large on proper scales 
up to 20Mpc at z = 5.5 and rapidly increases at earlier times 
before full reionization (Wyithe & Loeb 2006). The SDSS 
presently contains only a handful of quasars at a high enough 
redshift (z > 6.1) to probe a large path length of the z > 6 
IGM. From Fan et al. (2006a) and refs therein, the path 
length at z > 6 presently probed by the SDSS quasars is 
only ~ 40 proper Mpc (excluding the broad absorption line 
quasar SDSS J1048+4607, whose broad absorption lines may 
be confused with IGM absorption). Observations of SDSS 
quasars at z > 6.2 show evidence for significant cosmic vari- 
ance with a substantial range of foreground IGM absorption 
(Oh & Furlanetto 2005). More quasars are required to over- 
come cosmic variance at z > 6. 

Observations of high-redshift quasars are important not 
only for studying reionization, but also for the formation of 
supermassive black holes and their host galaxies. The high 
luminosities and broad line widths of the most distant quasars 
require black holes with masses greater than 10 9 M Q (Fan et 
al. 2001; Willott et al. 2003). Forming such massive black 
holes within the first billion years of the universe provides 
a challenge to models of galaxy formation, black hole for- 
mation and black hole growth (Bromley et al. 2004; Yoo & 
Miralda-Escude 2004; Shapiro 2005; Volonteri & Rees 2005; 
King & Pringle 2006). 

It is well established that there is a tight relation between 
the mass of the black hole and the spheroid in local galaxies 
(Magorrian et al. 1998). However the evolution of this re- 
lation with redshift is not well known and attempts are now 
being made to determine the relation up to z = 3 because it is 
a crucial component of galaxy evolution models (Shields et 
al. 2003; Adelberger & Steidel 2005; Alexander et al. 2005; 
Peng et al. 2006; McLure et al. 2006). It is also important 
to push this relation out to z ~ 6 where the earliest known 
black holes exist. Studies of millimeter radiation show that at 
least 40% of z ~ 6 quasars have high far-infrared luminosi- 
ties indicative of high star formation rates of ~ lOOOM yr _1 
(Wang etal. 2007). SDSS Jl 148+5251 has resolved CO emis- 
sion making it the only z > 6 quasar with a dynamical mass 
measurement (Walter et al. 2004). Curiously, the mass in- 
ferred is an order of magnitude less than expected if there is 
no evolution in the black hole - spheroid relation. 

To date, only 10 quasars have been discovered at redshift 
greater than or equal to 6 (Fan et al. 2006a and refs therein; 
McGreer et al. 2006). To discover more very distant quasars 
which can be used as probes of reionization and the early 
growth of massive black holes and galaxies, we are carry- 
ing out the Canada-France High-z Quasar Survey (CFHQS). 
This project uses the 1 square degree imager MegaCam at the 
Canada-France-Hawaii Telescope. The survey plans to cover 
~ 900 square degrees in 4 optical filters g',r',i',z' to a limit- 
ing depth of z' ~ 22.5. This is more than 2 magnitudes fainter 
than the SDSS survey leading to the discovery of correspond- 
ingly lower luminosity quasars. Lower luminosity quasars are 
expected to have lower black hole masses which could be due 



to the black holes being younger than those in SDSS quasars. 
The survey is carried out in areas of the sky which already 
have MegaCam imaging in some filters as part of the CFHT 
Legacy Survey 11 or RCS2 Survey 12 . As discussed below, the 
survey is also extremely effective at finding ultracool brown 
dwarfs. In this paper, we present the first high-redshift quasars 
to be discovered in the CFHQS. 

In the following section we describe how the quasars 
were discovered from imaging and spectroscopic observa- 
tions. Sec. 3 presents the properties of the four quasars. Sec. 4 
deals with observations at other wavelengths including mil- 
limeter observations. In Sec. 5 we discuss constraints on 
the ionization state of hydrogen at high redshift from these 
quasars. Sec. 6 presents our conclusions. 

All optical and near-IR magnitudes in this paper are 
on the AB system. Cosmological parameters of Hq = 
70 km s _1 Mpc -1 , Q, M = 0.27 and A = 0.73 (Spergel et al. 
2007) are assumed throughout. The convention for spectral 
indices, a„, is that f v oc v~ a " . 

2. OBSERVATIONS 
2.1. Imaging 

The four quasars whose discovery we report here were 
found in the first two years of optical imaging of the CFHQS. 
The total area covered at this time is « 400 square degrees 
to typical limits of z' w 22.5,/' ~ 24.5. However, because 
follow-up observations are not complete, we expect this re- 
gion to contain more quasars than discovered so far and no 
conclusions on the space density of low luminosity quasars 
and shape of the luminosity function at z = 6 are available at 
this time. 

We only briefly describe the processing and photometry of 
the imaging observations in this paper. A much fuller descrip- 
tion is given in Delorme et al. (2007). Pre-processing of the 
MegaCam images is carried out at CFHT using the ELIXIR 
pipeline. This removes the instrumental effects from the im- 
ages. We then run our own algorithms to improve the astrom- 
etry and check the photometry. Finally, we stack the images 
(for those cases where there is more than one exposure at a 
given position) and register the images in the four different 
filters. Photometry is carried out using our own adaption of 
Sextractor (Bertin & Arnouts 1996) which uses dual image 
multiple psf-fitting (Delorme & Bertin in prep.) to optimize 
the signal-to-noise ratio (S/N) of point sources. 

Candidate quasars are initially identified on the optical im- 
ages as objects which have very high i'—zf colors. In Willott 
et al. (2005a) we described in detail our method for identi- 
fying high-z quasars and brown dwarfs using MegaCam op- 
tical plus near-IR imaging. As shown in that paper, pho- 
tometric noise cause many M stars (with intrinsic colors of 
0.5 < i'—zf < 1) to be scattered into the region of the dia- 
gram at i' —z' > 1 .5 where we would expect to find only L 
or T dwarfs and quasars. The huge number of these M stars 
would require a lot of telescope time for complete follow-up. 
Therefore we now limit our survey to objects observed to be 
redder than i'—zf > 1.7, with stronger priority for follow-up 
given to redder objects. In Willott et al. (2005a), we showed 
that changing the criterion from - z' > 1 .5 to i' —z' > 1-7 
only loses us quasars at z < 5.9, but reduces the number of M 
dwarf contaminants by > 80%. Similarly, imposing a cut at 
i'-z' > 1.8 or even i'—zf > 2 would have very little effect on 

1 1 http://www.cfht.hawaii.edu/Science/CFHTLS 

12 http://www.rcs2.org 
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FIG. 1. — Color-color diagram showing the selection region for high-z 
quasars (to the lower-right of the thin solid line). The thin dotted line shows 
the optical criterion of i'—z' > 1-7 which also finds many L and T dwarfs. The 
thick solid line is the median cloned quasar track for quasars at 5.5 < z < 6.7 
with a mark every Sz = 0.1 (from Willott et al. 2005a). The thick dashed 
line is the median brown dwarf track for types M, L and T. The four squares 
with error bars are the four new z > 6 quasars. Colors of the brown dwarfs 
discovered by this survey will be presented in Delorme et al. (2007). 



the number of z > 6 quasars discovered. All of the quasars 
presented in this paper were observed to have i' -z 1 > 2. 

To separate quasars from brown dwarfs we employ near- 
infrared /-band imaging (e.g. Fan et al. 2001). Quasars at 
z < 6.5 appear relatively blue in z'-J, whereas L and T dwarfs 
are very red. Near-IR imaging was carried out at several tele- 
scopes and the observations, processing and photometry are 
described in Delorme et al. (2007). 

Fig. 1 shows the i' —z' vs z'-J color-color diagram for 
the newly discovered quasars. Also shown are lines corre- 
sponding to the expected median colors of quasars at red- 
shifts 5.5 < z < 6.7 and the median track for brown dwarfs 
discovered by 2MASS and SDSS. The expected quasar color 
track is based on "cloning" of 180 SDSS quasars at z = 3 (see 
Willott et al. 2005a for full details). Three of the quasars 
lie very close to the median quasar track. The exception is 
CFHQS J2329-0301 which is much redder in i'-z' and bluer 
in z' - J than expected given its redshift. The reason for this 
will be discussed in Sec. 3. Positions and photometry for the 
four quasars are given in Table 1 . 

20" x 20" images centred on the quasars are shown in 
Fig. 2. All quasars are clearly detected at z'-band, but most 
have only faint detections at i' and J. CFHQS J2329-0301 is 
the only quasar not detected at i'. CFHQS J 1509- 1749 is for- 
tuitously located only 12" from a very bright (R = 12.6) star 
and is ideally suited to natural guide star adaptive optics ob- 
servations. Finding charts in the z'-band over a wider field are 
presented in the appendix. 

2.2. Spectroscopy 

Optical spectroscopy of candidate quasars was carried out 
using the Gemini Multi-Object Spectrograph (GMOS; Hook 
et al. 2004) at the Gemini South Telescope and the Marcario 
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FIG. 2. — Images in the i', z' and J filters centred on the four CFHQS 
quasars. Each image covers 20" X 20". The images are oriented with north 
up and east to the left. 



Low Resolution Spectrograph (LRS; Hill et al. 1998) at the 
Hobby-Ebberly Telescope. The details of these observations 
are given in Table 2. 

The GMOS spectroscopy was carried out using the nod- 
and-shuffie mode to minimize sky subtraction residuals at the 
red end of the spectra. Most of the reductions of the Gemini 
spectra were performed using tasks in the IR AF Gemini pack- 
age, which were specifically designed for GMOS. The first 
step is removal of the bias level using a bias frame constructed 
from many (> 10) bias observations (with the same binning) 
carried out during the same month. Then the sky was removed 
from the spectra by subtracting a spectrum of the sky obtained 
through exactly the same light path. GMOS consists of three 
separate chips, so the next step was to mosaic together the data 
on the three chips. Each object was observed in more than one 
exposure, so all exposures of each object were registered in 
the spectral direction using the location of night sky emission 
lines (in unsubtracted images) and in the spatial direction us- 
ing the location of the quasar spectrum or that of a bright star 
located along the slit. The registered exposures were com- 
bined using bad pixel masks (including charge traps) and a 
sigma clipping algorithm. The exposures of CFHQS J0033- 
0125 and CFHQS J2329-0301 each had only four spectra to 
combine, so a few cosmic rays survived and were manually 
edited out of the final two-dimensional spectra. Wavelength 
calibration was achieved using the night sky spectrum from 
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Table 1 

Positions and photometry for the CFHQS quasars. 



Quasar 



RA and DEC (J2000.0) 



mag 



z mag 



/ mag 



CFHQS J0O331 1-012524 00:33:11.40-01:25:24.9 24.82±0.15 22.41±0.08 21 .58 ±0.20 2.41 ±0.17 0.74±0.22 

CFHQS J150941-174926 15:09:41.78-17:49:26.8 23.11 ±0.05 20.26±0.02 19.68±0.10 2.86±0.05 0.58±0.10 

CFHQS J164121+375520 16:41:21.64 +37:55:20.5 23.69±0.20 21.31 ±0.04 21.24±0.14 2.38± 0.20 0.07 ±0.15 

CFHQS J232908-030158 23:29:08.28 -03:01:58.8 > 25.08 a 21 .76 ± 0.05 21.56 ±0.25 > 3.32 0.18±0.25 



"Where not detected at > 2a significance, a 2a lower limit is given. 



Table 2 

Optical spectroscopy observations of CFHQS quasars. 



Quasar 


Redshift Telescope 


Date 


Resolving Slit Width Exp. Time 


Seeing 


^1450 




z 






Power 


(Arcsec) 


(s) 


(Arcsec) 




CFHQS J0033-0125 


6.13 


Gemini 


2006 Nov 26 


1300 


1.0 


3600 


0.6 


-25.03 


CFHQS J 1509- 1749 


6.12 


Gemini 


2006 May 26 


2100 


0.5 


7200 


0.7 


-26.98 


CFHQS J 164 1+3755 


6.04 


HET 


2006 July 26 + Sep 19 


400 


2.0 


5700 


1.9 


-25.48 


CFHQS J2329-0301 


6.43 


Gemini 


2006 Nov 26 


1300 


1.0 


3600 


0.7 


-25.23 



the unsubtracted images. The quasar spectra were then ex- 
tracted from the two-dimensional spectra. 

The LRS data reduction followed similar procedures. Cos- 
mic rays were removed from the individual frames prior to 
co-addition with L. A. Cosmic (van Dokkum2001). The wave- 
length scale was set and the distortion removed from the 
frames with a combination of cadmium and neon lines. A 4th 
order fit to the wavelength scale yields an accuracy of 1.0 A 
or better from 6500 -9 100 A. Incomplete fringing removal, 
beating with the night sky lines, causes some residual positive 
features in the final spectrum. These are described further in 
Sec. 3. 

For the GMOS and LRS spectra, a relative flux calibra- 
tion was applied using observations of a spectrophotometric 
standard star. Due to slit losses, absolute flux calibration was 
achieved by passing the spectra through the CFHT MegaCam 
z' filter + CCD response curves and normalizing to match the 
observed z' magnitudes. 

The optical spectra of the four quasars are shown in Fig. 3. 
In all four cases there is an unambiguous spectral break close 
to the Lya emission line, indicating the object is a z > 6 
quasar. All of the quasars except CFHQS J0033-0125 have 
strong Ly a emission lines. 

Note that the spectra have a wide range of S/N due to the 
spread of 2 magnitudes in brightness and different telescopes 
and exposure times used. Therefore we leave a detailed dis- 
cussion of the spectra to individual sections about each quasar 
in Sec. 3. 

Near-infrared spectroscopy of CFHQS J 1509- 1749 was 
performed using the Gemini Near Infra-Red Spectrograph 
(GNIRS; Elias et al. 2006) at Gemini-South on 2006 March 
14. Conditions were photometric with 0.8 arcsec seeing. 
GNIRS was used in the cross-dispersed mode with the 32 
lines mm -1 grating and the short camera. With a 0.675 arc- 
sec wide slit, this setup gives the full near-IR spectrum from 
0.8 to 2.4/im at a constant resolution of R = 900. The total 
integration time of 3300 s was split into 1 1 equal exposures to 
enable background subtraction and cosmic ray rejection. 

The GNIRS spectrum was reduced with the IRAF Gem- 
ini package. The individual frames were flat-fielded and then 



background subtracted using neighbouring frames with the 
target shifted along the slit by 3 arcsec. The spectra were 
rectified and then individual orders were extracted from the 
frames. The 1 1 individual frames for each order were then 
registered and combined using a bad pixel mask and cosmic 
ray rejection. Each order was wavelength calibrated with the 
argon arc lamps. For orders 3 to 6, the wavelength calibra- 
tion was refined using sky emission lines. The spectra were 
flux calibrated using a spectrum of an AO type standard star. 
One-dimensional spectra were extracted from each order and 
combined. Absolute flux calibration was achieved by passing 
the spectrum through the J filter response curve and normal- 
izing to match the observed J magnitude. The near-IR spec- 
trum of CFHQS J 1509- 1749 is shown in Fig. 4 and discussed 
in Sec. 3. 

3. NOTES ON INDIVIDUAL QUASARS 

3.1. CFHQS J2329-0301 

The spectrum of this quasar shows a very high equivalent 
width Ly a emission line. The line has a strong narrow peak 
at A = 9029 ±2 A. The FWHM of this narrow peak is only 
600 km s" 1 . Taking this narrow peak to be the systemic red- 
shift of the quasar (for reasons discussed below) would give a 
redshift of z = 6.427 ± 0.002. This makes CFHQS J2329-0301 
the most distant known quasar, surpassing SDSS Jl 148+5251 
at z = 6.419 ± 0.001 (Bertoldi et al. 2003b; Walter etal. 2003). 
The emission line has a broad base with a precipitous drop-off 
on the blue side at a Ly a redshift of z = 6.37. The spectral ex- 
tent of Ly a emission and transmitted flux in the Lyman forest 
will be discussed further in Sec. 5. 

We fit a power-law continuum curve by eye to the optical 
spectrum redward of Lya (> 9300 A) by fixing the curve to 
the observed /-band flux-density. Due to the low continuum 
level in this spectrum and relatively poor efficiency of GMOS 
at these wavelengths, the constraint on the power-law con- 
tinuum slope is fairly poor with spectral indices in the range 
0.5 < a„ < 1.5 acceptable. Planned infrared spectroscopy 
should yield a better value for a u and in this paper we adopt 
the 'maximum continuum' case where a u = 0.5 as is typical 
for unreddened quasars. 
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FIG. 3. — Optical spectra of the four newly discovered quasars. The expected locations of Ly a, Ly /3 and Ly 7 are marked with dashed lines. The two quasars 
with the highest S/N spectra have inset panels which show the Ly a and Ly f3 forests on an expanded scale. Note the huge variation in the strength and shape of the 
Ly a lines of the four quasars. All spectra are binned in 10A pixels. The spectrum of CFHQS J1641+3755 is contaminated by night sky residuals at A > 7600 A. 
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With the continuum level defined as above, we measure the 
rest-frame equivalent width of the observed Ly a line to be 
160 A. The main uncertainty on this measurement is that of 
the continuum level. Since we have adopted the maximum 
continuum level, this gives a lower limit to the equivalent 
width. We also note that because of absorption of the Ly a 
line, the intrinsic equivalent width before absorption will be 
even higher. A rest-frame equivalent width of > 160 A for 
Lya is very high compared with that measured for a com- 
posite spectrum of SDSS quasars of 93 A (Vanden Berk et 
al. 2001). Another way of comparing the Lya line to that 
of other quasars is the height of the peak of the line com- 
pared to the continuum. For CFHQS J2329-0301, the peak 
is 20 times the continuum level. This compares with a typ- 
ical peak height of only 3 times the continuum level in the 
SDSS composite. All 19 quasars found at z > 5.7 in the SDSS 
have peaks less than 10 times the continuum level (Fan et al. 
2006b). Therefore it appears that CFHQS J2329-0301 has an 
unusually strong narrow component to its Ly a line. The high 
Ly a EW explains the unusual location of this quasar in Fig. 1 
because Ly a falls in the z' filter. 

If the peak of the Ly a line suffered from significant absorp- 
tion, then the intrinsic peak would be even higher. We con- 
sider this unlikely given how high the Ly a peak and equiva- 
lent width are. Therefore we are confident that the measured 
Ly a redshift of z = 6 ATI ± 0.002 is close to the systemic red- 
shift of the quasar and adopt z = 6.43 ± 0.01 as the redshift and 
its uncertainty. We plan to obtain spectroscopy of the Mg II 
emission line to compare with the Ly a redshift, since Mg II 
usually shows only small offsets from Ly a (Richards et al. 
2002). We note that previous work on high-redshift quasars 
has usually defined the redshift via the Ly a peak (Fan et al. 
2006a and refs therein) and that the few high redshift quasars 
with Ly a and Mg II redshifts show no systematic offset be- 
tween them (Wyithe et al. 2005). 

3.2. CFHQS J0033-0125 

This is the faintest of the four quasars with z' = 22.41 which 
is close to the CFHQS magnitude limit of z' = 22.5. The 
spectrum therefore has low S/N. In contrast to CFHQS J2329- 
0301, CFHQS J0033-0125 has a very weak Ly a emission line 
which barely rises above the continuum. There is a narrow 
peak in Ly a at A = 8674 A which we tentatively identify as 
the systemic redshift giving z = 6.13, but the relatively low 
height of this peak above the continuum makes this redshift 
estimate quite uncertain. Plausible values for the redshift lie 
in the range 6.10 < z < 6.18. 

The equivalent width of the Lya emission is difficult to 
measure due to the uncertain continuum but it appears to be < 
10 A. With the exception of this narrow peak, CFHQS J0033- 
0125 resembles the lineless quasars discovered in the SDSS 
(Fan et al. 2006a and refs therein). We suspect that the low 
equivalent width of Ly a in this quasar is due to an unusually 
high number of dense pockets of neutral hydrogen close to the 
quasar. 

Despite the low S/N, it is apparent that there is a higher 
transmitted flux in the Ly a forest at observed wavelengths 
7500 A - 8000 A than at 8000 A - 8500 A. This is ob- 
served in other quasars at similar redshift (see for example 
CFHQS J 1509- 1749 on the same figure) as expected based on 
the increase of the Ly a forest density with redshift and pro- 
vides confidence that this object really is a z > 6 quasar and 
not an unusual broad absorption line quasar such as have been 
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FIG. 4.— Near-IR spectrum of CFHQS J1509-1749 obtained with GNIRS. 
Several broad emission lines are labeled. Also shown are two metal absorbing 
systems at z = 3.266 and 3.392. Note that this spectrum has been smoothed 
so the Mg II doublets, which are resolved in the original spectrum, appear as 
single lines. 



found in the SDSS (Hall et al. 2002). 

3.3. CFHQS J 1509- 1749 

This quasar is the brightest discovered so far in the CFHQS 
with z' = 20.26. This magnitude is similar to the faintest 
z ~ 6 SDSS quasars. The optical spectrum is the most sen- 
sitive of the four presented here resulting in the highest S/N. 
The spectrum shows a strong continuum break across the 
broad Lya emission line. This line has a peak at about 
8660 A (z = 6.12) but is highly asymmetric with consider- 
ably more flux on the red side of this peak than the blue 
side due to neutral hydrogen absorption. The Ly a emission 
line has a rest-frame equivalent width of 48 A which is typ- 
ical of high-redshift quasars (unlike CFHQS J0033-0125 and 
CFHQS J2329-0301). A more detailed discussion of this op- 
tical spectrum is given in Sec. 5. 

The near-infrared spectrum of CFHQS J 1509- 1749 is 
shown in Fig. 4. The gaps at 1.4^m and 1.9yum are due to 
opaque wavelength regions of the Earth's atmosphere. Emis- 
sion lines of Ly a, C IV, C III] and Mg II are visible along with 
some weaker lines. Also evident are strong absorption lines 
in the /-band, which are discussed in more detail below. 

The Mg II emission line is a low-ionization line that shows 
significantly less asymmetry and blueshift (relative to the 
narrow [O III] line) than high-ionization lines such as C IV 
(Richards et al. 2002). Therefore it provides the best in- 
dicator of the systemic redshift of a quasar at redshifts too 
high for any narrow lines to be observable. We have per- 
formed a gaussian fit to the Mgll line profile assuming an 
underlying power-law continuum. Although it is known that 
the continuum around Mg II is not flat due to many blended 
weak Fe II lines, this does not affect the determination of the 
peak wavelength of a fitted gaussian. The peak of the gaus- 
sian for Mgll is at 1. 9923 ±0. 0015 ^m giving a redshift of 
z = 6.118 ± 0.006. We will hereafter refer to the redshift as 
z = 6.12. Note that this redshift agrees with the peak of Lya 
in both the GMOS and GNIRS spectra. The C IV line is sub- 
stantially offset from the Mg II and Ly a peaks. The mea- 
sured C IV peak is at 1. 0982 ±0.0010,um which is a blueshift 
from Mg II of 1200 km s -1 . Blueshifts of this size are common 
among lower redshift quasars (Richards et al. 2002). 

3.3.1. Absorption lines 
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Despite the moderate S/N ratio (~ 10) of the near-infrared 
spectrum of CFHQS J 1509- 1749, two strong absorption sys- 
tems of low ionization are clearly detected at z = 3.266 ±0.001 
and 3.392 ±0.001. The spectral resolution is R ~ 900, thus 
the two Mg II doublets are well resolved. They also are fully 
saturated as both lines of each doublet have roughly the same 
equivalent width: w obs (2796,2803) = 5.0 and 21.0 A for the 
Z= 3.266 and 3.392 systems, respectively. The 3cr equivalent 
width detection limit in the /-band is w bs,iim = 4.0 A, thus 
even the weaker Mg II system is well detected, at about the 
4cr level. 

The z = 3.392 system is extremely strong, w res t(2796) = 
4.8 A, with associated absorptions from the Mg I singlet, 
w rest (2852) = 1.5 A, and Fell lines (2344, 2382, 2586, 
2600 A), w rest (2600) = 2.9 A. This set of absorption lines is 
characteristic of damped Ly a absorbers (DLAs): about 80% 
of the Mg II systems with w rest (2796) > 1 .0 A are DLAs (Rao 
et al. 2006) and, moreover, all systems with w rest (2852) > 
0.7 A are DLAs (Rao & Turnshek 2000). This system has 
all the characteristics of a DLA. The z = 3.266 system has an 
Mg II doublet with w lest (2796) = 1 .2 A, but no other absorp- 
tion lines are detected. 

The Mg II doublet is detectable in this spectrum over the 
redshift range 2.10 < z < 3.83. This yields an incidence of 
w res t(2796) > 1.0 A systems of dN/dz ~ 1.2 for this single 
line of sight. This value of dN/dz is three times greater than 
the results from Mg II absorption studies at z ~ 2 which give 
dN/dz = 0.40 ± 0.05 for w rest (2796) > 1 A systems (Nestor et 
al. 2005; Prochter et al. 2006). Given the high Poisson noise, 
we cannot say that there is a statistically significant excess 
in absorption towards CFHQS J1509-1749. We do not ex- 
pect significant cosmic evolution of dN/dz for strong Mg II 
absorbers from z = 2 to z = 3.4, because only a weak evolu- 
tion, if any, is observed at 1 .0 < z < 2.3. Therefore the appar- 
ent excess towards this quasar, if significant, must be due to a 
different cause, perhaps related to the higher quasar redshift. 

To determine if this discrepancy is statistically significant, 
one can use the published results of Keck-ESI spectra of 
z ~ 6 SDSS quasars. Strong Mg II doublets are detected to- 
wards SDSS J1044-0125 at z = 2.2786 and SDSS J1306+0356 
at z = 2.20 and 2.53, while there is no Mgll absorp- 
tion towards SDSS J0836+0054, SDSS J1030+0524 and 
SDSS Jl 148+5251 (Becker et al. 2001; Djorgovski et al. 
2001; Fan et al. 2003; White et al. 2003; Ryan- Weber et 
al. 2006). The two Mg II systems of highest z have clear, 
associated Fe II (2586, 2600) strong absorption, as could be 
seen from the quasar spectra, although this is not mentioned 
by the authors. No equivalent width is reported for any of 
the Mg II lines, but w b s (2796) should be of order 5 A or 
higher for the three systems. The sum of the redshift ranges 
for each quasar over which strong Mg II absorption is de- 
tectable equals Az = 2.79 which yields dN/dz ~ 1.08. A 
similar value is obtained when adding the two systems to- 
wards CFHQS J 1509- 1749 to the sample, dN/dz ~ 1.11 (or 
0.61 <dN/dz < 1.60 at the la level). Thus we see tentative 
evidence for an enhancement in the absorber incidence for the 
most distant quasars. 

One possible explanation for this is that some of the most 
luminous z ~ 6 quasars are gravitationally lensed (Wyithe & 
Loeb 2002; Comerford et al. 2002). The high incidence of 
Mg II absorption would then be related to above average mass 
density along the line of sight. However, searches for multi- 
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FIG. 5.— Rest-frame UV spectrum of CFHQS J1509-1749 plotted on a 
log-log scale. At A res t > 2000 A the spectrum is well fit by an unreddened 
quasar template (blue dashed line). But at A rest < 2000 A the spectrum is 
considerably reddened and can be fit by a quasar template reddened by SMC 
dust (red dotted line). Note that this SMC dust reddened quasar overpredicts 
the flux at A rcst > 2500 A. 



pie images and lensing galaxies in the highest redshift SDSS 
quasars do not show any evidence for lensing (Fan et al. 2003; 
Richards et al. 2004; Willott et al. 2005b). There is no evi- 
dence for multiple images from the ground-based imaging of 
the CFHQS quasars discovered so far (Fig. 2). 

3.3.2. Dust reddening 

The full near-IR coverage of the GNIRS spectrum allows 
one to investigate the continuum shape of this quasar. A 
power-law fit to the /-band continuum gives a spectral index 
of a v = 1 .00 ± 0.02. This is considerably redder than typical 
quasars which have a„ = 0.5 (Richards et al. 2003). How- 
ever, inspection of the H- and K-b&nd spectra show that the 
spectrum flattens here. 

Because of the complex shape of quasar spectra at A res t > 
2300 A due to blended Fe II lines and Balmer continuum (the 
"small blue bump"), we do not fit a power-law to these data 
but instead overplot a composite quasar spectrum constructed 
from spectra of SDSS quasars (Fig. 5). This shows that a typ- 
ical quasar spectrum fits the continuum well in the H- and 
K-\mnds but not at shorter wavelengths (A res t < 2000 A). An 
analysis of extinction curves of reddened SDSS quasars (Hop- 
kins et al. 2004) has shown that quasar dust is most often 
similar to the dust in the Small Magellanic Cloud (SMC). 
This extinction curve steepens significantly in the UV com- 
pared to that of typical Milky Way (MW) dust. We also plot 
on Fig. 5 the composite quasar spectrum reddened by enough 
SMC type dust (A v = 0.10) to match the /-band spectrum. Al- 
though this reddened quasar template fits well at < 2000 A, it 
predicts a steeper slope at > 2000 A than is observed. We 
observe essentially no reddening of the intrinsic quasar slope 
over the range 2000 A - 3200 A. We have checked the possi- 
bility of reddening by dust with known extinction laws (e.g. 
MW, LMC, SMC) in the z ~ 3 absorbers and found fits no 
better than SMC dust at the quasar redshift. 

One possible explanation for the sharp change in extinction 
at A res t ~ 2000 A is that the dust reddening this quasar has 
been created by type II supernovae. In our galaxy, the dom- 
inant dust production mechanism is in low-mass AGB stars 
(Gerhz 1989). But at z > 6 there is insufficient time (< 1 Gyr) 
for such stars to have evolved off the main sequence. There- 
fore core-collapse supernovae (whose progenitors have short 
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main sequence lifetimes) may be the main dust producers in 
the early universe. The heavy metal content of supernovae ex- 
plosions is expected to produce vast quantities of dust with a 
quite different composition and grain size distribution to that 
from AGB stars (Todini & Ferrara 2001). Recent observations 
of a z > 6 broad absorption line quasar and a z > 6 gamma ray 
burst have both shown signs of extinction curves which are 
similar to that inferred here (Maiolino et al. 2004; Stratta et 
al. 2007). However, without rest-frame optical observations, 
it is impossible to say if this is due to supernovae dust (which 
would have considerable, but constant, extinction at 2000 A 
- 3000 A; Maiolino et al. 2004) or dust with only very small 
grains (so no extinction at all at > 2000 A). 

We conclude that CFHQS J 1509- 1749 shows unusual dust 
reddening. Infrared observations probing the rest-frame op- 
tical are necessary to determine if this dust matches the pre- 
dicted extinction curve of dust created by type II supernovae. 

3.4. CFHQS J 1641 +3755 

The optical spectrum of this quasar (Fig. 3) shows a broad 
Ly a line with a sharp drop in flux on the blue side. The Ly a 
line has typical asymmetry with the peak close to the sharp 
drop. We assume the peak at A = 8555 A defines the systemic 
redshift of z = 6.04. The redshift uncertainty is ±0.04. The 
spectrum of CFHQS J 164 1+3755 is affected by night sky sub- 
traction residuals caused by incomplete removal of fringing 
which lead to excess positive flux. This explains the apparent 
high transmission of the Ly a forest for this quasar. Because 
of this we cannot with certainty define the continuum level of 
this quasar redward of Ly a. Given the fact that the quasar has 
az'-J colour typical of quasars (Fig. 1) we assume a spectral 
index of a v = 0.5 to calculate the absolute magnitude from the 
/-band magnitude. 

4. MULTI-WAVELENGTH DATA 

We have searched the NASA/IPAC Extragalactic Database 
(NED) at the locations of the four new quasars and found no 
coincident objects from any other survey. Two of the quasars, 
CFHQS J0033-0125 and CFHQS J 164 1+3755 are located in 
regions of the FIRST radio survey (Becker et al. 1995). We 
have checked the FIRST image cutouts and neither quasar 
is detected to a flux limit of w 0.5 mJy meaning that neither 
quasar is radio-loud. 

4.1. 1. 2 mm MAMBO observations. Dust and ultra-luminous 
IR starbursts. 

Follow-up observations of large samples of 2<z<6 
quasars at mm and sub-mm wavelengths have shown that a 
large fraction (30%) of optically luminous quasars are hyper- 
luminous infrared (L F i R > 10 13 L Q ) sources. In these sources, 
the far-IR luminosity is mainly related to the warm (40-60 K) 
dust, with estimated dust masses of few 1O 8 M . The heating 
of the warm dust appears to be dominated by the starburst ac- 
tivity of the quasar host galaxy and the implied star formation 
rates are of - lOOOMoyr" 1 (e.g. Omont et al. 2001; 2003). 
In a growing number of cases, warm and dense molecular gas 
is detected via CO emission lines and, in a few cases, in other 
species, revealing the presence of large reservoirs of molec- 
ular gas, the fuel of the star forming activity (see reviews by 
Cox et al. 2005 and Solomon & Vanden Bout 2005). The 
presence of such huge starbursts in phases of strong accre- 
tion of the quasars proves the simultaneity of major phases 
of growth of the most massive (elliptical) galaxies and their 



Table 3 

mambo 1.2mm photometry of cfhqs quasars. 



Quasar 


Redshift S L2 (mJy) 


CFHQS J0033-0125 


6.13 


1.13±0.36 


CFHQS J 1509- 1749 


6.12 


1.00±0.46 


CFHQS J1641+3755 


6.04 


0.08±0.46 


CFHQS J2329-0301 


6.43 


0.01±0.50 



super-massive black holes, and is an important clue for ex- 
plaining the black hole - spheroid relation. 

The detection of mm or sub-mm dust emission has been re- 
ported in eight z > 5.7 SDSS quasars (but only three at z > 6) 
with typical flux densities at 1.2 mm of a few mJy (Bertoldi 
et al. 2003a; Priddey et al. 2003; Robson et al. 2004; Beelen 
et al. 2006; Wang et al. 2007). Their properties and detec- 
tion rate, as well as the CO detection in J 1 148+525 1 (Bertoldi 
et al. 2003b; Walter et al. 2003, 2004), are comparable to 
other high-z strong starburst quasars. This shows the pres- 
ence of giant starbursts and vast amounts of molecular gas, 
with heavy elements, less than one billion years after the Big 
Bang. Jl 148+5251 is also the first high-z source where [C II] 
line emission was detected (Maiolino et al. 2005). 

It is important to extend the exploration of the millime- 
ter properties and star formation to a larger sample of z^6 
quasars, especially among the more numerous, fainter quasars 
in the CFHQS. A crucial question for building up the black 
hole - spheroid relation is whether the 1.2 mm emission and 
hence the star formation rate are strongly dependent on the 
UV luminosity, or only weakly as at 2 < z < 4 (Beelen 2004; 
Cox et al. 2005; Wang et al. 2007). Therefore we have ob- 
served at 1 .2 mm the four z > 6 CFHQS quasars. 

The observations were performed within the pool observing 
sessions at the IRAM 30m telescope in Winter 2006/2007, 
using the 116 element version of the Max Planck Millimeter 
Bolometer (MAMBO) array (Kreysa et al. 1998) operating at 
a wavelength of 1.2 mm (250 GHz). We used the standard on- 
off photometry observing mode, chopping between the target 
and sky at 2 Hz, and nodding the telescope every 10 or 20 
s (see e.g. Wang et al. 2007). On average, the noise of the 
channel used for point-source observations was about 35-40 
f" 1 / 2 mJy beam" 1 (where t is the integration time in seconds). 
This allowed us to achieve rms < 0.5 mJy for each of the four 
sources, with about 1.5 hr of telescope time per source. The 
data were reduced with standard procedures to minimize the 
sky noise with the MOPSIC package developed by R. Zylka. 

The results are shown in Table 3. There is no strong source 
with a flux significantly larger than 1 mJy. CFHQS J0033- 
0125 is detected at S/N> 3 with flux 1.13±0.36mJy. In ad- 
dition, CFHQS J 1509- 1749 has a very marginal (S/N ~ 2) de- 
tection. The mean flux in the four sources is 0.56±0.23mJy. 
This is an order of magnitude weaker than J 1148+5251. 
It corresponds nevertheless to a far-infrared luminosity 
> 1O 12 L0 and star formation rate of ~ 300M Q yr _1 , typical 
of ultraluminous infrared galaxies. This mean 1.2 mm flux 
for our four CFHQS sources is significantly smaller than the 
mean value that one can deduce from Tables 1 and 2 of Wang 
et al. (2007) for 17 z-6 SDSS quasars which is « 1.3 mJy. 
This difference could be related to the fainter UV luminos- 
ity of the CFHQS quasars. Observations of more CFHQS 
quasars in the future will allow tighter constraints on the 
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luminosity-dependence of far-infrared emission in the most 
distant quasars. 

5. IONIZATION STATE OF THE Z > 5 IGM 

The spectra of the quasars presented here are discovery 
spectra. They generally have modest integration times. Con- 
sidering the faintness of most of these quasars, detailed anal- 
ysis of the implications of these quasars for cosmic reioniza- 
tion is deferred until we have higher S/N and resolution spec- 
troscopy. In this section, we describe the information we can 
get from our observations on the ionization state of the high 
redshift universe. We only discuss the highest redshift quasar, 
CFHQS J2329-0301, and the brightest quasar, CFHQS J 1509- 
1749, since the other two quasars have very low S/N spectra. 
In addition, the systemic redshifts of CFHQS J0033-0125 and 
CFHQS J 164 1+3755 are poorly determined due to the lack of 
near-IR spectroscopy and low S/N of broad Lya. We note 
that poor sky subtraction can be a serious problem for studies 
of the high-redshift IGM. The spectra we use in this section 
are from GMOS in the nod-and-shuffle mode which is well- 
documented to provide excellent sky subtraction (Abraham et 
al 2004). 

5.1. Intrinsic quasar spectra 

In order to determine the neutral hydrogen absorption in 
these quasar spectra, we first need to determine the intrinsic 
spectra before absorption. This is not a trivial task given the 
large amount of absorption shortward of and around the Ly a 
emission line. As described in Sec. 3, we used optical/near-IR 
spectroscopy and/or photometry to determine the continuum 
power-law indices at A res , > 1250 A. CFHQS J 1509- 1749 has 
a v = 1 .0. The continuum for CFHQS J2329-0301 is not well- 
defined so we adopt a v = 0.5. Observations of a large sample 
of quasars at lower redshift show a continuum break in the re- 
gion 1200 A < A < 1300 A with a steeper power-law index of 
a v = 1 .7 at shorter wavelengths (Telfer et al. 2002). Therefore 
we adopt a v = 1.7 for both quasars at A rest < 1250 A. 

Because we are interested in absorption close to the quasar 
we also include broad emission lines of Ly (3, Ly a and N V. 
For CFHQS J 1509- 1749 we set the width of these lines in ve- 
locity to be equal to the measured C IV FWHM = 5600km s" 1 
(consistent with the Mg II FWHM). The peak height of the 
Ly a line is determined by eye and the Ly (3 and N V lines 
have the same fraction of the Ly a peak height as in the com- 
posite spectrum of Telfer et al. (2002). For CFHQS J2329- 
0301 the intrinsic width of the Ly a emission line is very un- 
certain. The spectrum shown in Fig. 3 has a strong narrow 
Ly a line on top of a broad base. Whether this observed pro- 
file is due to absorption of a single broad line or a combina- 
tion of broad+narrow lines is uncertain. We do not presently 
have near-IR spectroscopy for the measurement of other broad 
line widths. We adopt a single gaussian line with a width of 
FWHM= 3400km s -1 and peak equal to the observed Lya 
peak. The emission lines only have a significant effect very 
close to the quasar redshift, so this uncertainty does not affect 
the results presented. 

5.2. IGM optical depth 

The observed quasar spectra are divided by their intrinsic 
spectra as defined in Sec. 5.1 to determine the transmission 
as a function of wavelength. The transmission is then used 
to calculate the effective optical depth. When binning over 
multiple pixels, each pixel is weighted by the inverse variance 



Table 4 

igm absorption in bins of width az = 0. 15 

TOWARDS THE z = 6.12 QUASAR CFHQS J1509-1749 



Redshift 


Transition 


Transmission 




5.225 


Ly a 


0.1 10 ±0.004 


2.20 


5.375 


Ly a 


0.084 ±0.005 


2.48 


5.525 


Ly a 


0.051 ±0.005 


2.97 


5.675 


Ly a 


0.026 ±0.005 


3.64 


5.825 a 


Lya 


0.000 ±0.005 


>4.64 


5.825 b 


Ly/3 


0.097 ±0.022 


5.25 



"Optical depth is a 2a lower limit. 

b Transmission is after correction for foreground Ly a as described in the 
text and the optical depth is the equivalent Ly a optical depth using r^ ya = 
2-25 T Ly g. 



to reduce the effects of noise spikes due to bright sky emission 
lines. For CFHQS J2329-0301 it is evident that the low S/N 
of this discovery spectrum prevents us from placing strong 
constraints on the optical depth of the IGM along this line- 
of-sight. We detect flux in the region of Ly a absorption at 
z < 5.75, but not at higher redshifts. In a large bin covering 
the redshift range 5.35 < z < 5.75 we measure a Ly a optical 
depth of r = 2.67 ±0.17. In the higher redshift bin of 5.75 < 
z < 6.35 we obtain only a 2a limit on the optical depth of 
t > 3 .93 . This lower limit for the optical depth is considerably 
lower than the detections and limits for z > 6.2 SDSS quasars, 
so we do not discuss it further in this paper and await higher 
S/N spectroscopy. 

The spectrum of CFHQS J 1509- 1749 does have high 
enough S/N to investigate the optical depth of the IGM. We 
bin the spectrum in intervals of Az = 0.15 to get the effec- 
tive optical depth along this line-of-sight. We ensured that the 
highest redshift bin does not include flux that may be asso- 
ciated with the quasar proximity zone discussed in Sec. 5.3. 
Values of Ly a transmission and effective optical depth are 
given in Table 4. The uncertainties reflect only the photon 
noise, which is the dominant cause of measurement uncer- 
tainty. However, in most cases, this uncertainty is lower than 
the intrinsic cosmic variance (Fan et al. 2006b). 

The Ly (3 transmission is measured only for the highest red- 
shift bin. Lower redshift bins of Ly (3 would overlap with the 
higher order Lyman transitions. The Ly (3 region at z = 5.8 is 
contaminated by foreground Ly a absorption at z = 4.7. This 
foreground Lya absorption is corrected for statistically us- 
ing the curve of effective Ly a optical depth against redshift 
measured for a large number of lines-of-sight at z < 5.5. Fan 
et al. (2006b) quote this as r LyQ = 0.85[(1 +z)/5] 43 . Apply- 
ing this correction gives a value for the Ly (3 transmission of 
0.097 ± 0.022. The quoted uncertainty does not include the 
uncertainty in the foreground Ly a transmission. For this bin 
we detect zero flux in Ly a but there is a > 4(7 detection of 
flux in Ly (3. This is because the oscillator strength for the 
Ly (3 transition is considerably lower than for Ly a. 

Fig. 6 plots transmission against redshift for the five Lya 
bins and one Ly (3 bin. Also plotted are Ly a transmissions 
for SDSS quasars in Fan et al. (2006b). The solid curve 
is an extrapolation of the fit to z < 5.5 quasars described 
above. This curve is consistent with the increasing absorption 
at higher redshifts being solely due to the increase in number 
density of Ly a forest clouds for typical model IGM param- 
eters (Miralda-Escude et al. 2000; see also the model curves 



10 



Willott et al. 




5.4 5.6 5.8 
Redshift z 




FIG. 6. — Lyman transmission in the spectrum of CFHQS J1509-1749. 
Solid (red) squares show the Ly a transmission and the open (blue) square is 
the Ly f) transmission in bins of width Az = 0. 15. Small circles show the Ly a 
transmission from the sample of 19 SDSS quasars of Fan et al. (2006b). The 
solid line is the fit to the Ly a transmission of z < 5.5 quasars given by Fan et 
al. (2006b). The dashed line was obtained by Becker et al. (2007) by fitting 
flux PDFs of z < 5.4 quasars with a lognormal optical depth distribution. At 
z > 5.7 our data (and most of those in the SDSS) fall below the extrapolation 
of the Fan et al. curve. 



FIG. 7. — Effective Gunn-Peterson optical depth from the spectrum of 
CFHQS J1509-1749 (symbols as in Fig. 6). The Ly J3 point (open blue square) 
is the equivalent Ly a optical depth. The solid and dashed lines are the ex- 
trapolations from lower redshift quasars by Fan et al. (2006b) and Becker et 
al. (2007) (as in Fig. 6). The dotted and dot-dashed lines are the effective 
optical depth evolution in the simulations of Paschos & Norman (2005) and 
Gnedin & Fan (2006), respectively. 



in Fig. 7)). The Ly a data for CFHQS J1509-1749 lie close to 
the curve at z < 5.4 but fall well below it at higher redshifts, 
indicating a sharp evolution in the hydrogen neutral fraction 
compared to the expected extrapolation. 

Becker et al. (2007) have cast doubt on this form for the 
expected evolution in transmission. They showed that the 
Lya transmitted flux probability distributions of quasars at 
2 < z < 5 are better fit by lognormal optical depth distribu- 
tions than the distributions from theoretical models such as 
Miralda-Escude et al. (2000). Because the transmitted flux 
for the highest redshift quasars primarily depends upon the 
low optical depth tail of the distribution in voids, using the 
correct form for the optical depth distribution is important. 
The dashed curve in Fig. 6 shows the expected transmission 
evolution based on fits to quasars at z < 5.4 by Becker et al. 
It is clear that this curve falls much more rapidly than that 
of Fan et al. (2006b) and is consistent with the z > 5.8 data 
for no evolution in the hydrogen neutral fraction. However, 
the Becker et al. curve does not do a good job of fitting the 
data at 5.4 < z < 5.7, where it predicts transmission consid- 
erably lower than most of the SDSS data points. It is beyond 
the scope of this paper to investigate the cause of this dis- 
crepancy, but we note that Becker et al. used linear fits to 
evolution parameters, so it is perhaps not surprising that this 
linear fit cannot accomodate the data at all redshifts. 

Fig. 7 shows the optical depth evolution. For the highest 
redshift bin, the Ly a point is a 2a lower limit. For the Ly (3 
point, we plot the equivalent Lya optical depth assuming 
r LyQ = 2.25 T Ly/3 (Fan et al. 2006b; Gnedin & Fan 2006). The 
Ly a limit is consistent with the Ly (3 detection. The line-of- 
sight to CFHQS J 1509- 1749 shows one of the highest optical 
depth measurements at z ~ 5.8. As for the previous plot, it 
is clear that there is a sharp evolution away from the Fan et 
al. (2006b) low redshift curve at z > 5.4. Assuming standard 
IGM parameters (e.g. Fan et al. 2001), the optical depth evo- 



lution for CFHQS J 1509- 1749 is consistent with an increase 
in the IGM hydrogen neutral fraction of a factor of three from 
z = 5 to z = 6. 

Also plotted on Fig. 7 are two simulation results for the 
Gunn-Peterson optical depth evolution. The simulations of 
Paschos & Norman (2005; dotted curve) show no evolution 
away from the low redshift extrapolation in this redshift range 
because reionization is 99% complete by z = 6.4. Our data and 
that of the SDSS clearly show that the optical depth is greater 
than in this model at z > 5.7. 

The dot-dashed curve is the optical depth evolution for the 
L8N256 simulation of Gnedin & Fan (2006). In these simula- 
tions the overlap period of reionization occurs at z = 6.2. The 
emissivity at z = 6 ±0.1 in this model is set to match the SDSS 
observations, so we should expect agreement at this redshift. 
This simulation predicts a very sharp increase in the optical 
depth at z > 6. As is evident from Fig. 7, the actual evolu- 
tion for CFHQS J 1509- 1749 and SDSS quasars is somewhat 
more gradual than in this model. This suggests that the period 
of overlap is likely more extended than in the simple model. 
As highlighted by Gnedin & Fan (2006), the optical depth at 
the end of overlap is dominated by Lyman limit systems and 
simulations need to both resolve very small scales and cover 
very large scales. Both simulations gives r evolution during 
the epoch in which reionization is complete comparable to the 
Fan et al. (2006b) curve and much flatter than the Becker et 
al. (2007) curve. We defer a more detailed comparison of the 
observed optical depth evolution with models until we have 
high S/N spectra of more CFHQS quasars. 

5.3. Highly-ionized near-zones 

One of the difficulties of probing the hydrogen neutral frac- 
tion via quasar spectroscopy is that Lyman series absorption 
is extremely efficient. Because the density of the IGM in- 
creases at high redshifts, a neutral fraction (by volume) of 
only 2 x 10" 4 at z = 6 gives an optical depth at Ly a of r ~ 5 
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(e.g. Fan et al. 2006b). 

Therefore it is useful to identify methods which are more 
sensitive to high neutral fractions. One of these which has 
been extensively studied in recent years is using the size of 
the region around the quasar which it has ionized. This is 
variously referred to as an HII region, Stromgren sphere or 
highly-ionized near-zone. We will adopt the term near-zone 
as this does not imply that the surroundings are significantly 
neutral. 

For a quasar embedded in a uniform density IGM, the size 

-1/3 

of the ionization front from the quasar is proportional to / HI , 

N 1 / 3 and f^ 3 where /hi is the fraction of hydrogen in the IGM 

that is neutral, N is the number of ionizing photons emitted by 
the quasar per second and ?q is the age of the quasar activity 
(Cen & Haiman 2000; Wyithe & Loeb 2004). Note that the 
latter two terms could be combined into the total number of 
ionizing photons emitted so far by the quasar. Therefore a 
measurement of the size of the ionization front (from spec- 
troscopy) combined with knowledge of the number of ioniz- 
ing photons emitted would give the neutral fraction. 

This method has been employed to argue that the high- 
est redshift SDSS quasars are embedded in a highly neutral 
IGM (/hi £0.1; Wyithe & Loeb 2004; Wyithe et al. 2005; 
Mesinger & Haiman 2007). However there are many compli- 
cating factors such as the uncertain density distribution close 
to luminous quasars and pre-ionization by associated galaxies 
(Yu & Lu 2005). Uncertainty in the systemic redshifts can 
also be an issue leading to a typical uncertainty in the near- 
zone size of ss 1 Mpc for redshifts from the Ly a line. 

An additional uncertainty comes from relating observations 
to simulations. Defining the exact end of a near-zone from 
spectra is difficult. This is especially true for 5.7 < z < 6.2 
quasars which typically do not show Gunn-Peterson troughs 
immediately shortward of their near zones. In order to be 
able to define near-zone sizes for a wide range of spectra, 
Fan et al. (2006b) adopted the definition as the point in the 
spectra where the Ly a transmission first dropped to T < 0.1 
for spectra binned in 20A pixels. Most studies involving just 
those quasars at the upper end of the SDSS redshift distribu- 
tion (e.g. Wyithe et al. 2005; Yu & Lu 2005) define the near- 
zone size as corresponding to the onset of the Gunn-Peterson 
trough (typically T<0.01). Bolton & Haehnelt (2007a) dis- 
cuss many of the uncertainties associated with using near- 
zone sizes to estimate /hi. Chief among these is the fact that 
the observed size may not correspond to the actual size of the 
ionized region due to significant residual neutral hydrogen left 
behind the ionization front. 

As mentioned previously, the sizes of near-zones are ex- 
pected to be related to the quasar luminosity. For the classical 
case of an expanding H II region the near- zone size is propor- 
tional to N 1/3 . Bolton & Haehnelt (2007a) found the expected 
dependence to be N l l 2 for near-zones defined as T < 0. 1 . The 
CFHQS magnitude limit is 2.3 magnitudes fainter than that 
of the SDSS meaning that the CFHQS quasars will be on av- 
erage almost a factor of 10 lower in luminosity. Therefore 
the CFHQS near-zone sizes are expected to be between 2-3 
times smaller than those in the SDSS. Comparing the near- 
zones sizes for the final CFHQS sample to the final SDSS 
sample will allow a test of the luminosity dependence. 

In the upper and middle panels of Fig. 8 we show the trans- 
mission in the spectrum of CFHQS J 1509- 1749 close to the 
quasar redshift. The middle panel shows the Ly a and Ly (3 
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FIG. 8. — Lyman transmission close to the quasar redshift for 
CFHQS J1509-1749 (upper two panels) and CFHQS J2329-0301 (lower 
panel). The horizontal scale shows the distance in physical coordinates from 
the quasar. The upper panel shows the Ly a and Ly /3 unbinned spectra for 
CFHQS J1509-1749. The Ly (3 transmission has not been corrected for fore- 
ground Lya absorption and is shown vertically offset by +0.6 for clarity. 
The middle panel shows the same data but binned in 20A bins with Ly /3 
now corrected for foreground Ly a absorption. The lower panel shows 20A 
binned Lya transmission for CFHQS J2329-0301. Dotted horizontal lines 
show T = 0. 1 . 



transmission binned in 20 A bins. It is seen that the Ly a trans- 
mission first drops below T = 0. 1 at a distance of 4.3 Mpc, but 
then rises up to T = 0.2 in the next bin and then falls once 
more below T = 0.1 at 6.4 Mpc. This behaviour of T > 0.1 
peaks beyond the distance where T first drops below 0. 1 is 
also seen in the near-zones of some SDSS quasars (Fan et 
al. 2006b) and in the simulated quasar line-of-sight of Bolton 
& Haehnelt (2007a) and is likely due to to the existence of 
low-density regions within the quasar near-zone. Statistically 
significant transmission below the T = 0.1 level is detected in 
bins at distances up to 12.7 Mpc. Given these observations, 
what should we conclude about the size of the highly-ionized 
near-zone surrounding this quasar? 

To further understand the properties of the near-zone we 
also show the binned Ly (3 transmission in the middle panel of 
Fig. 8. These data have been statistically corrected for fore- 
ground Ly a absorption. Excluding the first few Mpc, there 
is statistically significant Ly (3 transmission all the way out to 
12.7 Mpc. In fact, the Ly (3 optical depth in this region is con- 
siderably lower than expected given the observed Ly a optical 
depth, which suggests the foreground Ly a forest at z ~ 4.9 
along this line-of-sight is unusually transparent. 

The upper panel shows the Ly a and Ly (3 transmission (this 
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time without statistical correction of the foreground Ly a for- 
est) of the full resolution spectrum without binning. There is a 
correspondence between positive features in Ly a and Ly at 
most places, although the correspondence in width and height 
of these features vary due to the varying foreground Ly a ab- 
sorption. The transmission feature at 12Mpc is clearly visible 
in Ly 0. From the weakness of this feature at Ly a we con- 
clude that the Ly feature is due both to a low-density region 
near the quasar and an unusually transparent patch at z w 4.9. 
We note, as emphasized by Songaila (2004), that if flux is de- 
tected at all at Ly a, it is better to adopt this measurement than 
to trust the apparently higher S/N Ly measurement because 
of uncertainty in the foreground subtraction. Nevertheless, 
the detection of positive flux in all Ly a bins out to 12.7 Mpc 
suggests that this is the true size of the zone impacted by the 
quasar ionizing flux. This is significantly greater than the 
6.4 Mpc that would result from imposing a T < 0.1 criterion. 

Throughout the above discussion we have assumed that the 
flux detected is not due to transparency caused by a very low 
neutral fraction (/hi < 10" 4 ) in the IGM beyond the quasar 
ionization front. There are a couple of facts supporting this 
assumption. First, there is significant flux detected in every 
Ly bin from 3 to 12.7 Mpc and positive (though not always 
significant) Ly a flux in every bin. Then at larger distances, 
the IGM appears very opaque in Ly a and Ly until there is 
some Ly transmission at ~ 20 Mpc. The transmitted flux in 
the region 6.4 to 12.7 Mpc corresponds to an effective Lya 
optical depth of r = 3.0. This is to be compared with the op- 
tical depth at similar redshift (this region is 5.92 < z < 6.02) 
and bin size along the lines-of-sight to all four highest redshift 
SDSS quasars of r > 4.3 (Fan et al. 2006b). If the transmis- 
sion at 6.4 to 12.7 Mpc is beyond the quasar ionization front, 
then the IGM here has an unusually high ionizing background 
over a substantial redshift range. 

The lower panel of Fig. 8 shows the Ly a transmission close 
to CFHQS J2329-0301. The spectrum does not have a high 
enough S/N to provide useful Ly data. The Ly a near-zone 
is remarkably similar to that of CFHQS J 1509- 1749 with the 
transmission first dropping below T = 0.1 at a distance of 
3.7 Mpc followed by a spike with T > 0.1 and then falling be- 
low T = 0. 1 again at 6.3 Mpc. Due to the lower S/N of the data 
we would not be able to detect any peaks at larger distances 
of comparable strength to those in CFHQS J1509-1749. 

The CFHQS quasars are mostly of lower luminosity than 
those of the SDSS so we must scale our near-zone sizes to a 
common luminosity before comparison with SDSS quasars 
and simulations carried out to match SDSS quasars. We 
scale the sizes to the equivalent for M1450 = -27 assuming 
R oc N 1 ? 3 . Because CFHQS J1509-1749 has M 1450 = -26.98, 
no luminosity-scaling is necessary. For CFHQS J2329-0301 
with M1450 = -25.23, the near-zone size of 6.3 Mpc goes up to 
10.8 Mpc after luminosity-scaling. 

We can compare the luminosity-scaled sizes of the CFHQS 
near-zones with those from the more luminous SDSS quasars. 
Bolton & Haehnelt (2007b) give the luminosity-scaled sizes 
of the Ly a near-zones measured from the four SDSS quasars 
at z > 6.1 (excluding the broad absorption line quasar 
SDSS J1048+4607) using a similar definition to us of the last 
pixel where T > 0.1. The four SDSS near-zone sizes are 2.6, 
4.8, 6.2 and 13.6 Mpc. Our two CFHQS sizes of 6.4 and 
10.8 Mpc are within the range of the SDSS quasars, some- 
what towards the upper end of the distribution. 

Although there are many uncertainties associated with us- 



ing the near-zone sizes to estimate the neutral fraction of the 
IGM that the quasar ionization front expands into, the finding 
of relatively large near-zone sizes in the two CFHQS quasars 
does allow us to determine some constraints. Firstly, as shown 
by Bolton & Haehnelt (2007a), there is little dependence upon 
the quasar lifetime if /hi ^ 0.1 and/or at ?q > 10 7 yrs because 
in these cases the observed near-zone size depends not on the 
expanding ionization front, but on the residual neutral hydro- 
gen left behind the ionization front. Both short lifetimes and 
high neutral fractions lead to small near-zone sizes, so we can 
safely assume that quasar lifetime is not a factor. Bolton & 
Haehnelt (2007b) simulated many quasar lines-of-sight with 
parameters of M1450 = -27 and fQ = 10 7 yrs. They found that 
for /hi > 0.3, the near-zone sizes are all < 4.5 Mpc. At 
/hi = 0.1, the average near-zone size is 4.5 Mpc with a range 
up to 6.5 Mpc. The only simulated near-zone sizes greater 
than 10 Mpc come from simulations with / H i < 10~ 3 . 

The luminosity-scaled sizes for the two CFHQS quasar 
near-zones are 6.4 and 10.8 Mpc. Therefore we can be quite 
confident that neither quasar is located in an IGM with /hi > 
0.3. For CFHQS J2329-0301 at z = 6.43 with a luminosity- 
scaled near-zone size of 10.8 Mpc, a very highly ionized IGM 
is inferred. Even if the quasar redshift has been overestimated 
by Sz = 0.02, then the observed near-zone size of 5.2 Mpc 
would still give a luminosity-scaled size of 9.0 Mpc. For both 
quasars, the quasar radiation is emitted into a substantially 
pre-ionized IGM. This indicates that these locations had al- 
ready been reionized to / H i < 0.3 at their respective redshifts 
of z = 6.12 and z = 6.43. Several authors have suggested that 
because luminous quasars are expected to form in rare over- 
dense regions, the surrounding IGM had already been pre- 
ionized by galaxies (Yu & Lu 2005; Alvarez & Abel 2007; 
Lidz et al. 2007). The possibility that quasar near-zones do 
not represent typical regions of the IGM must always be kept 
in mind when interpreting results such as these. However, we 
note that the lower luminosities of most CFHQS quasars, such 
as CFHQS J2329-0301, mean they are expected to be powered 
by less massive black holes hosted in less rare, lower mass 
dark matter halos, reducing the likelihood of pre-ionization 
by cluster galaxies. 

6. CONCLUSIONS 

We have presented the discovery of four new CFHQS 
quasars at z > 6. These were discovered from incomplete 
follow-up of less than half the final survey area and we expect 
many more quasars to be found in the next few years. We 
have discussed the observed properties of the quasars includ- 
ing emission lines, intervening absorption lines, continuum 
slopes and dust reddening. 

Sensitive millimeter continuum observations have been car- 
ried out with MAMBO at the IRAM 30m telescope. One of 
the four quasars has been securely detected. This strong dust 
continuum is likely due to a high star formation rate in its 
host galaxy. The mean 1.2 mm flux of the CFHQS quasars 
is substantially lower than the mean 1.2 mm flux of SDSS 
quasars at similar redshifts. This is indicative of a luminosity- 
dependence of the dust continuum emission. 

We have investigated the constraints that the discovery 
spectra can put on the ionization state of the IGM. For all but 
one quasar, higher S/N spectroscopy are required to perform 
a full analysis of the IGM properties. 

CFHQS J 1 509- 1 749 at z = 6 . 1 2 shows significant IGM opti- 
cal depth evolution from z = 5 to z = 6. The highest redshift bin 
at 5.75 < z < 5.90 has an effective optical depth of r Lya « 5, 
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which is one of the highest effective optical depths found at 
this redshift and consistent with a factor of three increase in 
/hi fromz = 5. 

Both CFHQS J 1509- 1749 and CFHQS J2329-0301 have 
large highly-ionized near-zones of at least 6 proper Mpc. 
These large near-zones are inconsistent with quasar ionization 
fronts expanding into a highly neutral (/hi <L 0.3) 1GM (Bolton 
& Haehnelt 2007b). Our results are broadly consistent with 
the (mass-averaged) value of / H i <~ 0.04 at z ~ 6.4 determined 
by Fan et al. (2006b) from evolution of SDSS quasar near- 
zone sizes. Our strongest constraint is from CFHQS J2329- 
0301 at z = 6.43 which has a luminosity-scaled near-zone size 
of 10.8 Mpc which is only consistent with / H i < 10" 3 . In ad- 
dition, the lower luminosity of this quasar compared to SDSS 
quasars makes it less likely that its surroundings were pre- 
ionized by cluster galaxies. 

Taking together the results of (i) a significant increase in 
the IGM optical depth at z « 5.8 (for CFHQS J 1509- 1749 and 
also observed in the SDSS) and (ii) the low neutral fractions 
inferred from the near-zones of CFHQS (and some SDSS) z > 
6.1 quasars, suggest that reionization is well under way before 
z = 6.4 but is still not complete by z = 5.7. This is a more 
extended period than results from most simple models where 
overlap occurs rapidly. 

Finally we note that cosmic variance is still a significant 
issue for the highest redshift quasars and that many more 
quasars at z > 6.1 are required to accurately determine the 
IGM evolution at z > 6. 
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APPENDIX 
FINDING CHARTS 

Fig. A9 presents 3' x 3' finding charts for the CFHQS quasars. All images are centred on the quasars and have the same 
orientation on the sky. These are the MegaCam z'-band images in which the quasars were first identified. MegaCam has gaps 
between the CCDs and these data were not dithered so the gaps remain and are evident in the upper two images where the quasars 
lies clos to the edge of a CCD. For all quasars except CFHQS J 1509- 1749, the Megacam data are a single exposure leading to 
many cosmic rays in the final images. 
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FIG. A9. — z'-band finding charts for the CFHQS quasars. 



